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Abstract 

The radii of nucleoli distribution and neutron skin in nuclei beyond the /3-stability 
line are studied within the extended Thomas-Fermi approximation. We show that 
the growth of neutron skin in unstable nuclei does not obey the saturation condition 
because of the neutron coat. The neutron coat indicates the possibility of giant neutron 
halo which is growing with moving away from the beta stability line. We demonstrate 
the presence of strong shell oscillations in the charge radius Rc and the relation of 
Rc to the isospin shift of neutron-proton chemical potentials AA = A n — A p for nuclei 
beyond the beta-stability line at fixed value of mass number A. 
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1 Introduction 

In the vicinity of the beta-stability line, the average changes in binding energy E and nuclear 
radius R with nucleon content obey the saturation properties. The volume part E vo \ of 
binding energy and the nuclear volume itself are proportional to the particle number A with 
E vo \ = —by A and R = tqA 1 ^ 3 , where by > and r and are constants. Both values of by 
and tq depend, however, on the isotopic asymmetry parameter X = (N — Z)/ (N + Z). This 
dependence comes from the difference in saturation bulk density, po ~ r o~ 3 ' for nuclei with 
different values of X. The saturation density po becomes smaller beyond the beta-stability 
line for neutron-rich nuclei where more neutrons are pushed off to form the " neutron coat" . 
One can expect that the growth of neutron skin in neutron-rich nuclei violates the saturation 
property R ~ A 1 ^ 3 for the nuclear radius providing an existence of neutron halo (giant 
neutron halo) effect pQ. 

In this paper we study the deviation of neutron distribution from the saturation behavior 
in neutron-rich nuclei. We study the influence of spin-orbit and Coulomb forces on the 

neutron, y/ (r^) , and proton, J (r^), root mean square radii as well as the relation of the 
shift a/ (r%) — y (rl) to the surface symmetry energy. We study the problems related to 
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the nucleon redistribution within the surface region of the nucleus and, in particular, the 
neutron coat and the neutron excess for the nuclei beyond the beta stability line. 

We combine the extended Thomas- Fermi approximation (ETFA) and the direct varia- 
tional method assuming that the proton and neutron distributions are sharp enough, i.e., 
that the corresponding densities p p (r) and p„(r) fall from their bulk values to zero within 
a thin surface region. In our consideration, the thin-skinned densities p p (r) and p n ( r ) are 
generated by the profile functions which are eliminated by the requirement that the energy 
of the nucleus should be stationary with respect to variations of these profiles. 



2 Direct variational approach 

We will use the ETFA which is one of practical realization of general Hohenberg-Kohn 
theorem [2] on the unique functional relation between the ground state energy and the local 
density of particles for any fermion system. The total kinetic energy of the many-body 
fermion system is given by the semiclassical expression [31 H] as follows 

E kin {p n , p p } = E kin {p g , VpJ = J dr e kin [p n (r), p p (r)], (1) 

where e kin [p n ,p p ] = e kin>n [p n ] + e kill)P [p p ], and 



ekin, q [p q , Vp 9 ] - ^ 



5 K ' Hq 36 p q 3 



Here p q is the nucleon density with q = n for neutron and q = p for proton. 

We will follow the concept of effective nucleon-nucleon interaction using the Skyrme-type 
force. The total energy functional for charged nucleus is given by 

Etot{pq, VpJ = E kiD {p q , VpJ + E SK {p q , VpJ + E c {p p }, (2) 

where E SK {p q , V p q } is the potential energy of iViV-interaction 

#sk{p<z, VpJ = J dr e pot [p n (r),p p (r)], (3) 

e pot [p n (r ), p p (r)] is the potential energy density and is the Coulomb energy. The potential 
energy ([3]) also includes the energy of spin-orbit interaction. Considering the asymmetric 
nuclei with X = (N — Z)/A 1, we will introduce the isotopic particle densities, namely 
the total density p + = p n + p p and the neutron excess density p_ = p n — p p with p_ p + . 
We apply the direct variational method (5] and assume the density profile functions p+(r) 
and P-(r) to be a power of the Fermi function as 

P+(r) = pof(r), p_( r ) = ^/(r)-ipo^A. (4) 
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Here, /(r) = [1 + exp [(r — R)/a}}~ v , the values p an d Pi are related to the bulk density, R 
is the nuclear radius, a is the diffuseness parameter and A is the parameter of neutron skin 
(see below). The profile functions p+{r) and p-(r) have to obey the condition of the neutron 
and proton number conservation. For the ground state of nucleus, the unknown parameters 
Po,i, R, cl, A, 7] and the total energy E tot itself can be derived from the variational principle 

5(E - \ n N - \Z) = 0, (5) 

where the variation with respect to all possible small changes of po,i, R, a, A and r\ is 
assumed. The Lagrange multipliers X n and X p are the chemical potentials for neutrons 
and protons respectively, and both of them are fixed by the condition of particle number 
conservation. 

As mentioned above, the parameter A in profile functions of Eq. (jlj) is related to the 
neutron skin. It can be easily seen from the derivation of the rms radii of the neutron and 
proton density distributions 

\[ijff = J f drr 2 p q (r) j J dr p q (r) . (6) 
Using Eqs. and @ one obtains the size of the neutron skin as 

Note that the evaluation of the variational conditions leads to an additional dependence of 
the variational parameters po,i ; R, a , A and r\ on the external parameters A and X. The 
value of A disappears in symmetric nuclei at X = and depends slightly on the Skyrme force 
parametrization. In case of the SkM forces we have numerically calculated the dependence 
of A on X for A = 120 and fitted it by the following formula 

A(X) w 0.90 X + 1.47 X 2 . (8) 

The parameter A is also related to the number, N$, of neutrons in surface region of the 
nucleus ("neutron coat"). Substituting Eqs. @j into condition of the particle conservation 
and using the leptodermous expansion, we obtain for the neutron excess N — Z the following 
expression 

N - Z w N v + N s , (9) 

where 

47T n / . , a , . a 2 \ 9 / . , a , , a? \ p n A 

N v = —R 3 (l + 3k ( V )- + 6«i(77)^ J Pi, N s = AttR 2 ^1 + 2k q ( V )- + 2^(^—1 ^A 

and Ki(r]) are the generalized Fermi integrals derived in Ref. [5]. The first term Ny ~ R? on 
the right hand side of Eq. [9]) is due to redistribution of the neutron excess within the nuclear 
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Fig. 1: Ns for neutron-rich nuclei in the vicinity of Sn (Z=50) (solid line). The dotted and 
dashed lines show the influence of the spin-orbit and the Coulomb forces on the neutron 
coat Ns- The arrow shows the position of Sn nucleus on the /3-stability line. 

volume while the second one Ns ~ R 2 is the number of neutrons within neutron coat. In 
Fig. [TJwe have plotted Ns for neutron-rich nuclei in the vicinity of Sn nucleus (Z=50) (solid 
line). Two additional lines show the influence of the spin-orbit and the Coulomb forces on 
the neutron coat N s . 

In general, the change of the radius R of nucleon distribution with the nucleon number 
A is caused by two factors. There is a simple geometrical change R oc A 1 / 3 . An additional 
change can occur due to the polarization effect (the bulk density distortion) with moving 
away from the beta-stability line. In particular, the size of neutron skin is sensitive to the 
symmetry and the Coulomb energies. We expand the total energy E tot (p , X) /A around the 
saturation density p ,eq and the isotopic asymmetry parameter X* on the beta-stability line 
as 

E tot (po,X)/A = E tot (po, cv X*)/A + t^Mpo - Po, cq ) 2 + ^-{X - X*) 2 (p - p 0j e q ), (10) 

18 Po,eq Po,eq 

where Ka is the incompressibility of finite nucleus 

d 2 E tot (p ,X*)/A 

K A = 9p , cq 7To (11) 

°P0 A, Po =po, cq 

and Pa is the partial pressure related to the symmetry and the Coulomb energies 

d 

P A = Pi eq^— [6v,sym(po) + &S.sym(po) A' 1/3 - a C {p )A 2/3 } , (12) 

°P° P0=P0, eq 
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where &v,s y m(Po) an d &s,s y m(Po) are the volume and the surface symmetry coefficients and 
a c(Po) — 3e 2 /20 (47rp /3) 1//3 . As seen from Eq. ffTUl) . the deviation from the beta-stability 
line (X ^ X*) implies the change in the bulk density po- The corresponding change is 
dependent on the incompressibility Ka and the partial pressure Pa- For an arbitrary fixed 
value of X, the equilibrium density po,x is derived by the condition 

= 0. (13) 

A,po=Po,x 

Using Eqs. ([TO]) and (fT3l) . we obtain the expression for the shift of bulk density (polarization 
effect) in the neutron rich nuclei 

Po,x = Po, cq -9^(X-X*) 2 . (14) 

The equilibrium partial pressure Pa is positive and thereby p 0j x < Po, eq , see also Refs. 

We point out that in general the sign of the equilibrium partial pressure Pa depends on the 

Skyrme force parametrization and this fact can be used for the Skyrme force selection [8]. 



—E tot (p ,X)/A 
opo 



3 Radii of nucleon distributions and neutron skin 

As noted above, the bulk density po,x is smaller for neutron-rich nuclei, more neutrons should 
be pushed off to enrich the skin providing the polarization effect. Thus, the rms radius a/T^U 
of neutron distribution does not necessarily obey the saturation condition ■sj (r 2 ) A 1 / 3 . 
As a consequence, the nuclei with significant excess of neutrons exhibit neutron skin, i.e., 
they are characterized by larger radii for the neutron than for proton distributions. The 
neutron coat Ns (see also Fig. [TJ) indicates the possibility of giant neutron halo which is 
growing with moving away from the beta stability line [lj. In Fig. [2] we have plotted the rms 
radii of neutron distribution from Eq. §6§ as a function of A. The deviation of v/(rf) from 
the saturation behavior ~ A 1 / 3 , obtained for the spherically symmetric nuclei, demonstrates 
the appearance of giant neutron halo when approaching the drip line. To extract a simple 
geometrical change of the radii we have made calculations with a step neutron distribution 
Pn(r) = po,n©(?" —Rn)i where the radius of the neutron distribution has saturation behavior 
R n = r ^ n A l l 3 . The results of the calculations are shown in Fig. |2]by the solid lines. As one 
can see from Fig. [2] the solid lines are very close to the beta-stability line. The difference 
between the dash-dotted and solid lines gives the value of the polarization effect. 

An occurrence of the giant halo in Na isotopes is shown in Fig. [3j We can see from 
this figure that the ETFA results agree quite well with the experimental data from [9]. The 
sensitivity of a/ (t^) calculation to the choice of the Skyrme force can be also seen. 

The results for the charge radius J (rfy are shown in Fig. HI The rms radius \J~if^) 
of proton distribution indicates the non-monotonic behavior due to the shell effects. Such 
behavior of J (r 2 ) correlates with A-dependence of the Coulomb radius Rc- To derive Rc 



5 



A 
V 



/23 Na 

... i i 



50 




208 



Beta-stability line 
Neutron rms radii 
Neutron drip-line 
Neutron rms radii (step distr.) 



_L 



_L 



100 



150 

A 



200 



250 



300 



Fig. 2: The rms neutron radii beyond the beta-stability line for spherically symmetric nuclei 
23 Na, 40 Ca, 91 Zr and 208 Pb. 




Fig. 3: The rms radius of neutron distribution in Na isotopes. 
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Fig. 4: The rms radius of proton distribution in Na isotopes. 

we use E tot (po,X) on the beta-stability line and establish an important relation for the 
chemical potentials X q (q = n for neutron and q = p for proton) beyond the beta-stability 
line. Namely, for fixed A, we obtain 



dE 
dN 



dE 
dZ 



d(E/A) 



N 



dX 



AX(X) 

where b* ym (A) = & VjSym (po )e q) + &s,sym(Po, cq ) A' 1 ^, 



4[b* aym (A) + e* c (A)](X-X*), (15) 



e* c (A) = 0A5Ae 2 /R, 



c 



and 



A, 



\dNj z 



A„= — 



dE 
dZ 



(16) 



(17) 



N 



On the beta-stability line, one has from Eq. (fl5|) that AX(X)x=x* = 0, as it has to 
be from the definition of the beta-stability line. We point out that for finite nuclei, the 
condition AA = on the beta-stability line is not necessarily fulfilled exactly because of the 
discrete spectrum of single particle levels for both the neutrons and the protons near Fermi 
surface. In agreement with Eq. (j!5p . the slopes of straight lines AA(X) allow us to derive 
the quantity 6* ym (A) + e* c {A). From the beta-stability condition AA(X) = one can also 
derive the asymmetry parameter X*(A). Using the definition of the beta-stability line as 



dE tot (p ,X)/A 



dX 



0. 



A, X=X* 



we obtain the symmetry energy coefficient b* ym (A) and Coulomb energy parameter e* c {A). 
Finally, using Eq. ffTB"]) . we obtain the A-dependence of the Coulomb radius Rc(A). 
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Fig. 5: The A-dependence of "experimental" Coulomb radius extracted from the experimen- 
tal data [ID] . 



The quantity d(E/A)/dX in Eq. ([15]) can be evaluated within the accuracy ~ 1/A 2 
using the finite differences which are based on the experimental values of the binding energy 
per nucleon B(N, Z) = —E(N, Z)/A. Namely, 

= $[B{N-1,Z + 1)-B{N + 1,Z-1)]. (18) 

A 4 

Since the difference (I18p is taken for AZ = —AN = 2, the pairing effects do not affect 
the resulting accuracy. Because of Eq. ( TIB]) , this procedure allows us to derive the "ex- 
perimental" value of the Coulomb radius Rc(A). The A-dependence of the "experimental" 
Coulomb radius Rq(A), extracted from the experimental data [TUJ, along the beta-stability 
line is presented in Fig. |3 The deviation of the Coulomb radius Rc(A) from the smooth 
A-dependence is mainly due to the shell oscillations. We point out that the shell oscillations 
in R* C (A) are correlated with the ones in A-dependence of the symmetry energy 

The size of the neutron skin ■sj (r^) — Wyf) i s illustrated in Fig. O The line has been 

obtained from Eq. ([7]), and the experimental data were taken from Ref. [9]. As seen from 

Fig. El the skin size \J (r%) — \/ ( r p) is primarily linear with the asymmetry parameter X. 



d(E/A) 
dX 



4 Summary 

We have applied the direct variational method within the extended Thomas-Fermi approx- 
imation with effective Skyrme-like forces to the description of the radii of nucleon distribu- 
tions. In our consideration, the thin-skinned nucleon densities p p (r) and p n (r) are generated 
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Fig. 6: Isovector shift of nuclear radius for Na isotopes. 

by the profile functions which are eliminated by the requirement for the energy of the nucleus 
to be stationary with respect to variations of these profiles. The advantage of the direct 
variational method is the possibility to derive the equation of state for finite nuclei: depen- 
dence of the binding energy per particle or the pressure on the bulk density p . We have 
evaluated the partial pressure Pa which includes the contributions from the symmetry and 
the Coulomb energies. The pressure Pa is positive driving off the neutrons in neutron-rich 
nuclei to the skin. 

Using the leptodermous properties of the nucleon densities p p (r) and p n ( r ), we have es- 
tablished the possibility of giant neutron halo in neutron-rich nuclei. The effect of giant halo 
increases with moving away from the beta stability line. In Fig. [2] this fact is demonstrated 
as deviation of the rms radii of neutron distribution from the saturation behavior ~ A 1 / 3 in 
the nuclei beyond the beta-stability line. 

The average behavior of the rms radii of nucleon distributions J (r^) and the size of 
neutron skin are satisfactorily described within the extended Thomas- Fermi approximation. 
The sensitivity of the calculations of J (r^j to force parametrization can be used for Skyrme 
force selection. The charge radius of proton distribution shows the strong shell oscillations 
with mass number. We have demonstrated the relation of the Coulomb radius to the isospin 
shift of the neutron-proton chemical potentials AA = A n — A p for nuclei beyond the beta- 
stability line at fixed value of mass number A. 
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